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Abstract 

First searches for the coherent dissociation of relativistic oxygen nuclei into four a particles 
are reported. It is shown that reactions of this type are characterized by a significantly lower 
decay temperature than the conventional multifragmentation of residual projectile nuclei. The 
momentum spectra and correlations of a panicles are not reproduced by the simple statistical 
model of direct fragmentation. The possibility that the oxygen nucleus undergoing fragmentation 
acquires a nonzero angular momentum in the collision process is discussed. 
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Multifragmentation of relativistic projectile nuclei provides unique insight into their struc- 
ture: as charged secondaries are detected without any thresholds, the corresponding tran- 
sitions can be observed at the lowest possible values of 4-momentum transfer. Extremely 
peripheral, coherent, inelastic nuclear collisions in which the target nucleus acts on the disso- 
ciating projectile nucleus as a discrete particle - that is, the former does not change charge, 
does not undergo breakup, and is not excited - are especially promising in this respect. 
Investigations of such collisions between high-energy particles were initiated largely by the 
milestone theoretical results by Pomeranchuk and Feinberg [T]. Note, however, that these 
results equally apply to collisions in which high-energy nuclei take the place of particles 
as projectiles, while the products of nuclear fragmentation appear as final states instead of 
elementary particles produced |2]. 

As in coherent particle production on nuclei, processes of two types - a diffractive tran- 
sition mediated by a Pomeron and a Coulomb transition mediated by a virtual photon - 
occur in the above nuclear reactions. The latter is expected to be dominant for colliding 
nuclei with high electric charges. Typically, inelastic coherent reactions have rather high 
energy thresholds. For example, the threshold incident momentum for coherent dissociation 
of a projectile nucleus A with mass Mq into n fragments with masses mi{i = 1, ...,n) is 
estimated as 

p^'"" ^ (Mo5^/V/i)A (1) 

where n is the pion mass, B is the mass number of the target nucleus that coherently interacts 
with the projectile A, and A = Yl^=i ^« ~ mass defect with respect to the 

dissociation channel considered. In particular, for the reaction ^^O — )■ 4a, the momentum 
threshold amounts to several hundred MeV/ c per nucleon, approaching the relativistic region 
for heavy targets. 

Only a few studies devoted to the coherent dissociation of relativistic projectile nuclei 
have been performed thus far. This paper reports the results of first searches for the coherent 
dissociation process 

160 ^ 4a (2) 

at po = 4.5 GeV/c per nucleon in nuclear emulsion (Em). Earlier, the dissociation of oxygen 
nuclei into four a particles was investigated at nonrelativistic energies (below 90 MeV per 
nucleon; see [3111] and references therein), which probably fall short of the coherent threshold 
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at least for the heavy target nuclei. 

Stacks of nuclear emulsion BR-2 irradiated by a beam of ^^O ions with a momentum of 
4.5 GeV/c per nucleon from the JINR synchrophasotron were used in searches for events 
satisfying criteria necessary to select reaction (2). 

To estimate the mean range A of ^^O nuclei in emulsion for reaction (2) without addi- 
tional charged secondaries, a slow-fast scan along the track was carried out in a part of 
emulsion sheets. Over the primary track of total length 375.2 m, we found 12 pure events 
of reaction (2). This corresponds to a A value of Sl.Sly^o^ m. Prom this, the cross section 
for the so-called average emulsion nucleus {{A) ^ 45) was estimated at 10 mb/nucleus. 

To maximize the statistics of events we are interested in, we made use of a special proce- 
dure based on a fast area scan in a large number of emulsion sheets. The scan was performed 
over bands normal to the incident oxygen beam. The distance between the bands was 0.5 cm. 
In the scan, we selected close groups comprising 3 to 4 gray tracks of similar ionization. Once 
such a group has been located, search along the group axis was carried out in the reverse 
direction. In comparison with scanning along the track, this strategy resulted in more than 
a tenfold increase in the scanning efficiency for reaction (2). 

In this way, we found and measured as many as 641 events that involve four well-identified 
relativistic fragments with z — 2 in the final state and which show evidence neither for the 
emergence of additional charged particles, nor for the excitation of the target nucleus. These 
events are further analyzed in the present study. 

Figure 1 shows the distribution in the square of the transverse-momentum transfer to 
the ^^O nucleus through interaction with an emulsion nucleus. This quantity can be used 
to estimate the square t of the 4-momentum transfer to the oxygen nucleus splitting into n 
fragments. This is done with the aid of the relation 

n 

\t'\ = |t-r^'^(^m,)| ^ ql (3) 

i=l 

Indeed, in the coherent diffractive production of the fragments (as in elastic diffractive 
scattering), the distribution in t' and in t as well, because the factor exp[t"""'(^"^-^ rrii)] that 
takes into account the minimum 4-momentum transfer squared t'^^"'{Y17=i'>TT'i) required for 
the production of n free fragments is insignificant here must have the simple exponential 
form 

da/dt' oc exp{-a\t'\), (4) 
4 
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FIG. 1: Distribution of ^^O 4a events in = —t'. A fit to the functional form (6) is shown 
by a curve. 

where the slope of the diffractive peak a is expressed in terms of the radii of the projectile 
and target nuclei, Ra and Rb, as a ^ {Ra + -Rb)^/4. On average, the resulting longitudinal 

component of the 3-momentum transfer q is much less than the transverse component. It 
follows that —t' ~ q^; hence, the distribution in the transverse momentum squared must 
have the Rayleigh form 

da/dq^ oc exp{—aq^) (5) 

with (qx) ~ h^^'^/{Ra + Rb)- We assumed that there are no neutral fragments in the 
events under study. In this case, we have = Ylt=i where pr is the a-particle 
transverse momentum. Its absolute value pt was estimated by the formula pt — 4:pQsin9, 
where po — 4.5 GeV/c and 9 is the polar emission angle. A straight hue corresponds to 
distribution (5) on the exponential scale used in Fig. 1. 

It can be seen from Fig. 1 that the experimental distribution in q^ is not consistent with 



the functional form (5) (x^/NDF = 3.2). However, it can be described by the sum of two 
Rayleigh distributions, 

da/dq^ — aexp{—aiq^) + (1 — a)exp{—a2qT), (6) 

with two different slope parameters Oi and 02- A maximum-likelihood fit of (6) to experi- 
mental data yields (xV^DF = 0.9) 

ai = 19 ± 2, 02 = 4.2 ± 0.4, a = 0.66 ± 0.06. (7) 

Once the nuclear composition of emulsion and experimental uncertainties in Qt, which effec- 
tively broaden the distribution in qt, are taken into account, the resulting values of Oi and 
02 are compatible with the assumption that the first and second terms in (6) correspond, 
respectively, to the diffractive dissociation ^^O — > Aa on an emulsion nucleus as a discrete 
particle and dissociation on a nucleon. At the same, any attempts at estimating the relative 
contributions of these two dissociation channels (that is, the slope parameters a) and their 
cross sections are seriously impeded by the fact that events with a recoil proton arc not 
observed in our experiment, by the possible contribution of Coulomb dissociation, by the 
area scan used, etc. 

Along with the complete sample of 641 pure events of the type ^^O — > 4q;, we will 
henceforth consider the subsample of 428 (~ 67% of the complete sample) events that 
satisfy the condition qt < 0.32 GeV/c, which is equivalent to t' < 0.1 (GeV/c)^. This 
value of t' (see Fig. 1) may be taken for an approximate boundary between the coherent and 
incoherent events in the selected sample. Provided that the coherent reaction (2) is indeed 
responsible for the first term in (6), this separation enables us to constrain from below the 
characteristics of the two reactions. In other words, we may hope to obtain some reliable 
estimates by comparing the characteristics of the complete sample and of the subsample 
with -t' < 0.1 (GeV/c)2. 

Figure 2 shows the transvcrsc-momcntum-squarcd distributions of a particles from reac- 
tion (2) for the complete sample and for the subsample with —t' < 0.1 (GeV/c)^. The 
corresponding root-mean-square values {pxY^"^ are 167 ± 4 and 145 ± 4 MeV/c. The two 
distributions display deviations from the Rayleigh form, showing high-energy tails that are 
due to particles with high pr- These deviations cannot be explained by the complex com- 
position of emulsion, in view of the well-known result from mathematical statistics that an 
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FIG. 2: Inclusive distribution of a particles for (a) the complete event sample and (6) events 
with —t' < 0.1 (GeV/c). The straight lines are the Rayleigh distributions with (p^) corresponding 
to the observed values. 

arbitrary superposition of any number of Rayleigh distributions with different a values is a 
Rayleigh distribution as well. 

It is also well known that the laboratory transverse momenta of fragments are effectively 
overestimated because of the transverse motion of the nucleus undergoing fragmentation (the 
so-called bounce-off effect discussed, for example, in [5]). In our data, this effect manifests 
itself in the distribution in the azimuthal angle Sij = arccos(pT.pTj /pt,PTj ) between the 
transverse momenta of two a particles produced in the same event of dissociation (2) (see 
Fig. 3). This distribution contradicts the phase-space prediction (the curve in the same 
figure). In particular, for the observed azimuthal asymmetry 

A = {N,^,>./2 - iV.,,<V2)/iVo<..,<. (8) 
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FIG. 3: Distribution in the aa azimuthal angle Eij for the complete sample of O — t- 4a events. 

of the distribution dcr/deij, we have -0.01 ± 0.02, while the value predicted by the law of 
transverse-momentum conservation is l/{Na - 1) = 0.33 [B]. This means that the fragmen- 
tation of the oxygen nucleus must be described in terms of the quantities defined in its rest 
frame. 

For an exclusive reaction like (2), the recipe for going over to the rest frame of the nucleus 
undergoing fragmentation is straightforward. (In the following, all the quantities that refer 
to this frame are labeled with asterisks.) Under the assumption that the projectile nucleus 
is scattered at a small angle, a high-precision approximation to the transverse momenta of 
a particles is given by the formula 

4 

i=l 

Figures 4 and 5 show, respectively, the p^. and e*j distributions of a particles for (a) the 
complete sample of events (2) and for (6) the subsample with t' < 0.1 (GeV/c). The 
main quantitative characteristics of these distributions are given in the table. These are the 
mean values (p^) and {p*^Y^'^-, the azimuthal asymmetry A* [see formula (8) for ejj], and 
the azimuthal collinearity B* defined as 

B* = (iV£*X7r/4 + ^4>37r/4 " N^/i<e*.<37r/4) /No<e*^<^. (10) 

The data presented in Figs. 4 and 5 and in the table lead to the following conclusions. 

(1) As might have been expected, the mean values (p^) are substantially smaller than 
(Pt) (by some 30% for the complete sample of ^^O — > 4a). It is noteworthy that the (p^) 
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FIG. 4: As in Fig. 2, but for the a transverse momentum defined in the ^^O rest frame. 

values for all events and for those with —t' < 0.1 (GeV/c)^ differ only shghtly. 

(2) As in the laboratory frame, the distributions show deviations from the Rayleigh 
form for both sets. It should be emphasized that within errors, the relative contribution of 
the high-p^ tail for the complete sample does not differ from that for the subsample with 

-t' < 0.1 (GcV/c)2 (23 ± 4 and 21 ± 4%, respectively). 

(3) The e*j distributions also disagree with the functional form 

da/ds* ~ —(1 + cicose* + C2Cos2e*), (11) 

which follows from the statistical picture of the direct nuclear decay into the observed a 
particles. Assuming that in the c.m.s., each component of the a 3-momentum obeys the 
normal distribution n{0, a) (this automatically leads to the Rayleigh form of the distribution 
da/dp^) and taking into account energy- momentum conservation in the decay, we find that 
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FIG. 5: Distributions in e*j for (a) the complete event sample and (6) events with 
—t' < 0.1 (GeV/c)^ The curves represent distribution (11) with coefficients (12). 

the coefficients Ci and C2 in (11) are related to the azimuthal asymmetry A* and the azimuthal 
coUinearity B* by the equations 



ci = -{n/2)A* = -(7r/2)(iV, - 1) 



-1 



C2 = (7r/2)S* = (87r/25)(7V, - 1)-^, (12) 

where (in our case) A^q, = 4 [6]. The predicted forms (11) with coefficients (12) are represented 
by the curves in Fig. 5, and the values of A* and B* are given in the table. The discrepancy 
between the predictions and experimental data is clearly seen. 

(4) In relation to the predictions based on transverse-momentum conservation, the ex- 
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TABLE I: 









Quantity 






Event sample 


MeV/c 


VeV/c 


A* 


B* 


a* 


13* 


Experiment, all events 


121 ± 2 


145 ± 3 




n 97 -1- n (19 

U.Z < It U.UZ 


-0.23 lb 0.01 


0.22 lb 0.02 


Experiment, 


115 ±2 


134 ±4 


n 97 -1- n n9 

— U.Z / It U.uZ 


n "in -A- n n9 

U.OU It U.UZ 


-0.23 lb 0.01 


0.25 lb 0.02 


events with t' < 0.1 (MeV/c)^ 














Statistical model, ^^0— > 4q 


120 




-U.oo 


n n7 

U.U ( 


-0.26 


0.06 


Statistical model. 


120 




-u.oo 


nil 

U.ii 


-0.27 


0.09 


16Q_^8gg_|_Q, _|_ Q, _^ 4q, 














Statistical model. 


120 




-u.oo 


n 1 n 

U. iU 


-0.27 


0.08 


leQ-j.SBe+'^Be-^ 4a 














Statistical model. 


120 




-0.33 


0.07 


-0.27 


0.06 
















Statistical model. 


120 




-0.33 


0.08 


-0.27 


0.06 


16q^12q* ^ ^ ^8Be+a + a ^ 4a 















perimental distribution da/de*j shows a higher degree of colhnearity between the a-particle 
transverse momenta in the transverse plane (see the observed excess of experimental 
values over theoretical results for Sij ~ and Eij ^ vr in Fig. 5; the values of B* are given 
in the table). 

Let us now discuss the results. We will use the Feshbach-Huang-Goldhaber statistical 
theory of prompt fragmentation [71 18] . In this theory (which is widely employed in analyses 
of fragmentation of excited nuclei), the measured values of (p^) can be related to the 
temperature of the oxygen nucleus decaying through the channel (2). In energy units, the 
temperature is estimated as 

KT = -^(a^/m^), (13) 

where m^r is the nucleon mass, and aj^ is the variance of the momentum distribution of 
intranuclear nucleons; the latter governs the variance of the momentum distribution for 
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arbitrary fragments through the so-called parabolic law 

4 = alAF{A-AF)/{A-l). (14) 

Here, A and Ap are the mass numbers of the nucleus undergoing fragmentation and the 
fragment under consideration, respectively. In our case, we have A = 16, Ap = A^ = 4, 
and 0"^ = {p^)/2. Using the experimental values of (p^) (see table) for the complete 
sample of events (2), we obtain kT = 3.7 MeV. For the events with t' < 0.1 (GeV/c)^ 
and with t' > 0.1 (GeV/c)^, we arrive at kT = 3.2 and 4.8 MeV, respectively. 

The above estimates indicate that the coherent dissociation ^^O — )■ 4a is characterized 
by a lower value of kT {kT 3 — 3.5 MeV) than conventional multifragmentation in 
reactions of the type A + B ^ a + X that are investigated in inclusive experiments 
(see, for example, [3 I5HTT]). The last statement remains valid even if we consider that in 
many experimental studies, the decay temperature was overestimated because of erroneously 
using the laboratory momentum characteristics of a particles without taking into account 
the transverse motion of the system undergoing fragmentation. The value that we obtained 
for kT is also significantly lower than the nucleon binding energy in the oxygen nucleus. On 
the other hand, such a low value of temperature corresponds to a small energy-momentum 
transfer to the nucleus undergoing fragmentation, which is typical of coherent processes. If 
nonstatistical tails in the distributions shown in Fig. 4 are due to the decays of massive 
intermediate states (the cascade decays of ^^O are discussed below), the true kT value is 
even lower than that quoted above. 

Let us consider this issue in more detail. First, we address the question of whether there 
is pt collinearity in individual events. For this purpose, we supplemented the analysis of 
the inclusive characteristics of azimuthal asymmetry and collinearity {A* and B*) with the 
calculation of the asymmetry a* and collinearity /3* in individual events by the formulas |6] 

4 

a* = {cose*^} = J2cose*j/No,{N^-l){i ^ j), (15) 

4 

a* = {cos2e%) = ^ cos24./iV„(iV, - 1)(2 ^ j), (16) 

(— (A^ — 1)"-*^ = —(1/3) < a*,/?* < 1). Assuming that the components of are normally 
distributed and that the transverse momentum is conserved in the decay, we find that the 
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FIG. 6: Distribution in the azimuthal collinearity /3* for an individual event. The solid and dashed 
curves represent the results of calculations made in the statistical theory for reaction channels (18) 
and (19). respectively. 

mean values (a*) and of random variables a* and P* are given by [S] 

(a*) = _(vr/4)(iV,-l)-i = -0.26, 

= (47r/25)(iV« - l)-2 = 0.056. (17) 

The distribution in the collinearity P* is plotted in Fig. 6 for the 428 ^^O — > Aa events 
with —t' < 0.1 (GeV/c)^. The observed mean values (a*) and {(3*) are presented in the 
last two columns of the table. 

To evaluate the contributions of various ^^O cascade decays to the observed tendency 
of the product a particles to diverge collinearly in the transverse plane, we performed a 
simulation of all possible cascade channels of dissociation according to the simplest statistical 
model [Bj. Thus, we considered the transitions 

^^O — )• 4a (the direct channel), (18) 

^^O ^ ^Be + a + a ^ 4a, (19) 
^^O ^Be + ^Be Aa, (20) 
^^O ^ + a -> 4a, (21) 
^^O a ^ + a + a ^ 4a. (22) 
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For each decay fragment from reactions (18)-(22), we assumed that each 3-momentum com- 
ponent defined in the rest frame of the decaying nucleus or of the unstable intermediates 
state obeys a normal distribution n(0, o"^), the variance o"^ being dependent on the fragment 
mass according to the parabolic law (14). After the transformation from the rest frames 
of intermediate states to the c.m.s. of the ^^O nucleus, we were left with the single model 
parameter a% [see (14)]. This parameter was determined by imposing the requirement that 
the mean transverse momentum {p"^) for the final-state a particles in each channel (18)- (22) 
be equal to its measured value ~ 120 MeV/c (see table). For each of the reactions (18)-(22) 
simulated according to the Monte Carlo method, the values of A*, B*, (a*), and are 
shown in the table. The computed distributions dN/d(5* are illustrated in Fig. 6 for chan- 
nels (18) (the direct reaction) and (19) (under the above assumptions, the latter leads to 
the maximum coUinearity of the vectors in the final state). 

The results shown in the table and in Fig. 6 demonstrate that the observed effect of 
collinearity cannot be attributed entirely to the cascade decays of ^^O into a particles. If 
we go beyond the model from [8j by assuming that the temperature of primary ^^O decays 
is substantially higher than that of the decays of intermediate fragments, a high degree 
of collinearity of the a-particle transverse momenta in the final state can be obtained. In 
this case, however, theoretical calculations fail to describe other experimental characteristics 
(mean values (p^), azimuthal asymmetries, etc.) under any assumptions about the relative 
contributions of channels (18) - (22). 

Therefore, the only plausible way to explain the high values of B* and is to assume 
that a nonzero angular momentum is transferred to the ^^O nucleus in the collision. Needless 
to say, a comprehensive analysis of conclusions that can be drawn on basis of this assumption 
is beyond the scope of the present study. We look forward to testing it with increased 
statistics that will be accumulated in the near future. 
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